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We investigate geodesies in specific Kundt type N (or conformally flat) solutions to 
■ Einstein's equations. Components of the curvature tensor in parallelly transported tetrads 

are then explicitly evaluated and analyzed. This elucidates some interesting global properties 
of the spacetimes, such as an inherent rotation of the wave-propagation direction, or the 
character of singularities. In particular, we demonstrate that the characteristic envelope 
qq ■ singularity of the rotated wave-fronts is a (non-scalar) curvature singularity, although all 

\Q ' scalar invariants of the Riemann tensor vanish there. 

o ' 
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O 

1 Introduction 

o 

In 1961 a general class of solutions to Einstein's equations which admit a non-expanding, shear- 
^jfj| free, and twist-free null congruence was introduced and described by Wolfgang Kundt [1-3] (for 

recent reviews see, e.g., [4], [5]). This contained some previously known families, in particular the 
pp -waves [6-9] or the Nariai [10] and Bertotti-Robinson [11-13] universes. Within his large class, 
Kundt also identified new solutions of various types, among them specific type ./V spacetimes 
representing plane-fronted gravitational waves which were geometrically different from famous 
pp -waves. These so-called Kundt waves are either vacuum solutions or contain pure radiation. 
Their generalization to a non-vanishing value of the cosmological constant A was later found by 
Ozsvath et al [14] (see also [15-18]). Very recently, these solutions have been further generalized 
in [19] where a complete family of Kundt waves of type III (which admit A and/or pure radiation) 
has been derived and classified. Geometry of the wave surfaces and their envelope singularity 
has also been analyzed. For particular choices of the metric functions one recovers more special 
type N Kundt waves and conformally flat non-expanding spacetimes. 

In fact, the family of conformally flat pure radiation Kundt metrics attracted attention 
few years ago [20-27]. In general, they contain no invariants or Killing or homothetic vectors, 
and thus provide an interesting exceptional case for invariant classification of exact solutions. 
It was demonstrated by Skea [25] that, to distinguish the Wils [20] metric within the more 
general Edgar-Ludwig solution [22] , it is necessary to go as far as but no further than the fourth 
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derivative of the curvature tensor. This result supports the conjecture that any space-time can 
be uniquely characterized using the derivatives of the curvature tensor in which no higher than 
the fourth derivative is ever required. 

Even more interestingly, it has been recently demonstrated [28,29] that for the above Kundt 
spacetimes all curvature invariants of all orders identically vanish. They may thus play an 
important role in string theory and quantum gravity since there are no quantum corrections to 
all perturbative orders [30]. 

Despite the fact that particular aspects of these spacetimes have been studied for decades, 
many questions concerning their geometrical and physical properties still remain open. For 
example, analysis of the geodesic deviation [31] revealed the structure of relative motions of test 
particles in the Kundt gravitational waves but it is not yet clear whether the geodesic motion 
itself does exhibit a chaotic behaviour, analogous to the chaos in pp -waves [32]. 

We consider here vacuum Kundt type N waves but some results also apply to conformally flat 
pure radiation solutions of this non-expanding class. The purpose of our paper is to contribute 
to understanding of their global structure, namely the character of singularities. For principal 
reasons this is not an easy task since all curvature invariants — which would unambiguously 
identify physical singularities — vanish. Thus, to analyse singularities in the Kundt spacetimes 
(summarized here in section it is necessary to study the behaviour of free test particles, which 
is done in section |21 In the subsequent section |1] we find frames parallelly transported along 
timelike and null geodesies, and in sectional we project the curvature tensor onto these tetrads. 
We employ these results for discussion of the singularities in the final section H3 

2 The Kundt spacetimes 

The metric of type N or conformally flat Kundt spacetimes, which are either vacuum or contain 
pure radiation, can be written in the form (see, e.g., [5]) 

ds 2 = 2d(dC-2Q 2 dudv + Fdu 2 , (1) 
Q = C + C, F = 2Q 2 v 2 -QH , 

where H (£, £, u) is a function of the spatial coordinates £, £, and of the retarded time u. A simple 
transformation v = wQ~ 2 puts the metric into the Kundt canonical form [1,4] 

ds 2 = 2dCdC-2dn (dw + W dC + WdQ + Hdu) , (2) 

with W = — 2itfQ~ 1 , 7i = —w 2 Q~ 2 + \Q H (the pp-w&ves corresponding to W = and 7i = H 
in (j2J will not be investigated here). In the natural null tetrad 

k = d v , 1 = {F/2)Q- A d v + Q' 2 d u , m = 9 C - , m = 8 ( , (3) 

(k is the quadruple principal null direction) the only non- vanishing Weyl and Ricci scalars are 

*4 = \ Q- 3 , $22 = \ Q~ 3 # C c • (4) 

The spacetimes are thus conformally flat when = 0, i.e. if H has a special form 

H = A + A 1 ( + A l C + A 2 C(, (5) 

where A4 are arbitrary complex functions of u. These spacetimes containing pure radiation 
characterized by $22 = ^2/(2Q 3 ) were considered by Kundt [1-3] and others [20-27]. For A2 = 
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Figure 1: The wave surfaces u = uq = tan(a/2) with x > are tangent half-planes to expanding 
cylinder (for T = const.; Left), cone (for Y = const.; Middle), or circle (for T,Y = const.; Right) 
in Minkowski background coordinates. At any time T each wave surface is bounded by the line 
where it touches the cylinder x = 0. For different values of a, the successive surfaces are rotated 
around the cylinder centered on the y-axis which expands with the speed of light (four surfaces 
u = uq are indicated; Left and Middle). As the circular caustic envelope x = expands (Right), 
the tangent half-planes propagate in the spatial direction perpendicular to each wave surface. 



the matter is absent, and Hq can be set zero by a suitable coordinate transformation [18]. Such 
solution with H = corresponds to flat Minkowski space, see transformation (|11|) below. 
Vacuum type N Kundt waves arise when $22 = 0, f 4 / 0, i.e. for 

H = f((,u) + f((,u) , (6) 

where / is an arbitrary function of £ and u, holomorphic in £, such that 7^ 0. These solutions 
(generalized also to admit a cosmological constant A) were presented in [14,15,18,19]. 

For our analysis of geodesies and singularities it is convenient to rewrite the Kundt solutions 
in real spatial coordinates x and y, related to the complex coordinate ( by 

( = (x + iy)/V2. (7) 

The metric Q obviously takes the form 

ds 2 = dx 2 + dy 2 - Ax 2 du dv + A(x 2 v 2 + xG) du 2 , (8) 

where G = -H/(2y/2). For vacuum spacetimes the function G(x, y, u) satisfies the Laplace 
equation, G :XX + G m = 0. The Christoffel symbols in coordinates (v, x, y, u) = (0, 1, 2, 3) are 

T° 01 = l/x, T° 03 = -2v, r° ls = -(G/x), x , r° 3 = -{G/x) tV , 

Tl 3 = 4v 3 + ivG/x-(G/x) tU , r| 3 = - W - 2G - 2xG, x , (9) 

^03 = , r| 3 = —2xG^ y , rf 3 = l/x , Tgg = 2v , 

and all independent non-vanishing components of the curvature tensor read 

-R1313 = — 2xG tXX , -R1323 = —2xG jXy , -R2323 = —2 x G, yy ■ (10) 

Let us also recall the geometry of the Kundt spacetimes in the weak field limit. It has been 
shown previously [26] that the transformation 

X = x(l + 2uv) , Y = y , 

Z = x [v — u(l + uv) ] , T = x [v + u(l + uv) ] , (11) 
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puts the Kundt metric (JSJ to 



xG(x,y,u) 
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ds 2 = ds 2 + - 2 ^ (1 + u 2 ) dT - 2udX - (1 - u 2 ) dZ , (12) 

where ds 2 , = — dT 2 + dX 2 + dY 2 + dZ 2 is the standard form of flat space corresponding to G = 0. 
In (|12j) . x, y and u should be expressed using the inverse of namely 



±V^ 2 + ^ 2 -T 2 , y = Y, 



X =f v 7 ^ 2 + Z 2 - T 2 T + Z 

U = , f = ± ; . (13) 

T + Z + Z 2 -T 2 V ' 

In the weak field limit (for small G) when the metric (|12j) is a perturbation of Minkowski 
background [26], each wave surface u = uq = const, corresponds to a hyperplane 

— T + sin a X + cos a Z = , y arbitrary , (14) 

where a = 2arctan«o- As ito increases from — oo to +oo then a goes from —it to +ir, and these 
wave surfaces form a family of null hyper planes which rolls all around the cylinder 

X 2 + Z 2 = T 2 , Y arbitrary . (15) 

For T > its radius is expanding at the speed of light. In fact, it follows from (|13j) that this 
caustic surface, formed as the envelope of all rotated wave surfaces, corresponds to a singularity 
x = in the metric ((SJ), see also [33]. 

To obtain a unique foliation (|14|) of the region outside the expanding cylinder (|15[) by the 
wave surfaces u = u$ G (— oo,+oo), it is necessary to restrict these to half-hyperplanes. This 
is achieved by considering only the range x > (or x < 0). In both cases the wave surfaces 
are localized outside the cylinder ()15|) given by x = which is formed as their envelope as 
indicated in figure ^ The restriction x > also resolves the ambiguity in the signs of (|13|) which 
corresponds to the invariance of the metric © under the reflections x — > —x, G — > —G. 

Natural questions arise concerning the character of the envelope x = of the rotating family 
of wave fronts, and the region inside it which is not covered by the Kundt coordinates. Is 
this envelope singularity only a coordinate one or is it a "physical" singularity? Is it possible 
to extend geodesies from the outside region across x = 0? And what tidal forces would the 
corresponding geodesic observers feel? We concentrate on these open problems in our paper. 

3 Geodesies in the Kundt spacetimes 

Considering the Christoffel symbols © the geodesic equations read 

x = -Axvii + (4xv 2 + 2G + 2xG, x )u 2 , (16) 
y = 2xG tV u 2 , (17) 
ii = -2(x/x)u- 2vii 2 , (18) 

(where the dot denotes differentiation with respect to an affine parameter r) plus a complicated 
equation for v. Instead of it, however, we take the normalisation condition of four-velocity: 

x 2 + y 2 - 4x 2 uv + Ax 2 [y 2 + G/x^j u 2 = e, (19) 
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Figure 2: Typical behaviour of timelike geodesies in the Kundt wave spacetime n = 2. The 
initial values are x = y = u = v = l, x = ii = — 1, y = 1, v = —2.75 at r = 0. 



where e denotes the character of the geodesic, namely 

{—1 ( timelike 

for < null geodesies. (20) 

+1 t spacelike 

Specific geodesic motion given by equations (|16|) - ()19J) depends on the particular spacetime 
function G(x, y, u). We concentrate here on type N vacuum Kundt waves given by the expres- 
sion ©. If the function f((,u) is linear in £, this generates only a flat space, cf. (|1J). The 
simplest radiative spacetimes thus arise when / = d{u)Q n , n = 2,3, 4, .... For a real function 
d(u), in coordinates x,y this corresponds to 

G (n) = c(u)Re{(x + iy) n } . (21) 

An arbitrary function c(u) determines the profile of gravitational wave. Various sandwich and 
impulsive waves can thus be constructed [34]. However, in our contribution we only consider 
spacetimes for which the function G is independent of u, i.e. c(u) = const. Without loss of gen- 
erality we assume c = 1. Indeed, employing the scaling freedom of the metric (|HI), u — » u = An, 
v — > v = A -1 v, G — > G = A~ 2 G, the factor c can be set to unity by an appropriate choice of A 
(G < corresponding to c = — 1 is related to G > with c = 1 by the reflection x — > —x). 

Even after these simplifications are made, it is difficult to find explicit solutions of equations 
(|16[) - (|19j) . To obtain more general insight into the geodesic motion it is necessary to perform 
numerical simulations. For example, typical behaviour of a test particle in the Kundt spacetime 
with G^ 

= x 2 — y 2 is indicated in figured It reaches the envelope singularity x = in a finite 
value of the proper time r, with the remaining coordinates decreasing. However, for different 
initial values the coordinate x may grow to infinity. Discussion of our numerical studies will 
be presented elsewhere [35]. Instead, we concentrate here on derivation of families of explicit 
geodesies which could be used for analytic investigation of the singularities. 

3.1 Explicit geodesies in the hypersurface y = 

The equations simplify if we restrict to yeodesics in the hypersurface y = 0. In fact, it is a 
natural assumption: in the background Minkowski coordinates this corresponds to Y = which 
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is a plane perpendicular to the symmetry axis of the cylindrical envelope singularity, see figure^ 
The expressions (|21|) for (with c = 1) and its derivatives on y = simplify to 



G (n) 



G^ 



nx 



n-l 



n{n) 
,2/ 



. 



(22) 



The equation (|T7j) is then identically satisfied by y(r) = 0. Remaining equations (fl6|) . (fl9|) . (fl8|) 

take the form 



x + 4:xvii - [4xv 2 + 2(n + l)x n ] u 2 



x 2 



n-l\ „\2 



4x 2 tt?) + 4x 2 (v 2 + x 

u + 2(x/x) u + 2t;u 2 



, 
. 



(23) 
(24) 
(25) 



From the relation (|25[) we easily express v (assuming and substitute into (|23j) and (|24[1. 

We obtain a system for the functions u and x only, 



3^ + 2(n + ^x™ -1 ?! 2 + 3- = , 
it x 



2- - 3^ + 4x n - 1 ii 2 
u it 



x x e 
— + 4 - = — 

X X X 



Subtracting 1)27(1 from (|26l) we derive an important relation 

e + 2(n - l)x 



2 x 



= XX + x 2 



' t+1 n 2 . 



(26) 
(27) 

(28) 



It is possible to express u and substitute into (|2^|) but the resulting fourth-order non-linear 
equation for x is very complicated, and additional assumptions are thus necessary. The term on 
the right-hand side of (|28[) indicates that the relation between u and powers of x is essential for 
the character of solution. For example, if x n+1 n 2 — > then it becomes (x 2 ) ~ 2e, which implies 



x/er 2 + Cir + C 



(29) 



where Co,Ci are constants. Such geodesies are possible: figure 01 shows null geodesies (e = 0) 
obtained by numerical simulations. Linear behaviour of the function x 2 near x = is obvious. 




1 1.5 
proper time X 




1 1.5 
proper time X 



Figure 3: Typical null geodesies approaching the singularity at x = in the Kundt spacetime 
n = 2. The initial values at r = are x = 0.5, 1, 1.5, and u = —0.45, —1.2, —2.17, respectively. 

To make the analysis more systematic let us assume a power-law relation between u and x, 

ii = Ax a , (30) 
where a,A/0 are constants. Putting this into (|28[) we get equation which can be integrated to 



e + 



AA 2 { 



n 



2a + n + 3 



1) x 2a+n+l + ^_ 



(31) 
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with C being the constant of integration. By substituting into (|26|) we derive the relation 



e Ti + AA 2 T 2 x 2a+n+1 + T 3 ^ = , 



(32) 



where the constant parameters are 

Ti = 
T 2 = 
T 3 = 



a{a + 2) , 
a + n + 1 r 



2a + n + 3 



a(n — 1) + 2n 



(33) 



1 



(2 + a) 



The equation (|3*2*|) is exactly satisfied if e = 0, C = 0, T 2 = 0, or if e = ±1, C = 0, 2a + n + 1 = 0, 
i.e. for three possible choices of the parameter a, 



(A) 
(B) 
(C) 



a 
a 
a 



-(n + 1) , 
-2n/(n - 1) 
-(n + l)/2 . 



(34) 
(35) 
(36) 



These cases enable three families of exact geodesies, denoted as A <n) , B <n) and C ( "', which will be 
discussed in detail in section l3~3l Before this, however, we first investigate a general asymptotic 
behaviour of solutions to (|31|1. (|3*2*|) — > or x — > co. 



3.2 Asymptotic behaviour of geodesies 

We start with asymptotic solutions as x — > 0. To avoid the divergence in the term x~ 2 in (|32|) 
it is necessary to set either C = or T3 = 0. All possibilities, summarized in tabled are: 

Asl. The region x — > is reached when C = 0, e = 0, and 2a + n + l = fc>0. These 
conditions imply x n+1 u 2 — > 0, i.e. near x = the solution approaches (|29|) . 

As2. For a = —2 we obtain T\ = 0. The region x — > can only be approached asymp- 
totically when C = 0, n > 3, e = 0, and the right-hand side of (|31|) is positive. This 
is a particular solution corresponding to the choice k = n — 3 in lAslI 

As3. The choice a = — 1 implies T3 = 0, and only for e = the equation (|32j) is fulfilled 
as x — ► 0. The right-hand side of equation (j31jl is then positive and it can be solved 
if C > 0. This asymptotic solution exactly corresponds to ()29|) . 

As4. We can also make T3 = by choosing a = —3. Remaining conditions are the same 
as in the case IAs3l but we have to assume n > 5 to reach x — > 0. In contrast to 
IAs31 T2 = 2 — n < 0. This solution also corresponds to ()29|) asymptotically. 





c 


a 


Eq. (jSU 


Eq. (JHU) 







|(fc-n-l) 


x fe ~ 0,fc> 









-2 


x n ~ :i ~ 0,n > 3 


i; 2 = \A l x n -' A 




> 


-1 


a-n-i ^ 


x 2 « Cx- 2 




> 


-3 


x n ~ 5 ~ 0, n > 5 


i 2 « Cx-^ 



Table 1: Asymptotic solutions as x — > 0, assuming u = Ax a . All these geodesies are null (e = 0). 
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Explicit forms of x(r) corresponding to the above four possibilities are summarized in tabled 
Notice that the solutions lAslI with < k < 2 reach x = in the finite value To of the affine 
parameter, whereas for k > 2 they approach this singularity as r — > oo. The solution |As2l exhibits 
an analogous behaviour, according to whether 3<n<5orn>5. 





explicit solutions as x — > 




x ~ 


|2-A;|yg^(r-ro) 


2 

2-fc 




2 

x ~ [ 5 — n\A (r — tq)] 5 ~ n 




x « 


2VC ; (t-to)" 


i 

2 




x sa 


'2VC(r-r )" 


1 

2 



Table 2: Explicit asymptotic forms of x(r) for null geodesies as x — > 0. 

We similarly investigate asymptotic geodesies in the region x — > oo. It is not now in general 
necessary to set C = since the term Cx -2 in l|31|). (j32j) becomes negligible for large values of 
x. Specific possibilities, summarised in table 01 are: 

As5. The region x — > oo is reached asymptotically if2a + n + l = A;<0 and e = 0. The 
case k = —2 is, however, forbidden. 

As6. The choice a = —2 implies T\ = 0, which requires n = 2 and e = to reach the area 
x — ► oo. This is a particular case of lAs5l for k = — 1. 

As7. If a = — 3 then T3 = 0. Asymptotic solutions x — > 00 now require n = 4 with e = 0. 
The alternative case n = 2 corresponds to IAs5l for A; = —3. 

As8. The choice (A) given by pi)) implies T 2 = 0. We thus have to put e = 0. 

As9. Choosing (B), i.e. Q35|) . we also obtain T 2 = which implies e = 0, 

AslO. Finally, the class (C) given by (|36|) (corresponding to A; = 0) is the only one which 
admits e = ±1. The Eq. (|32|) requires 2^4 2 = — e(n — 3)/(n 2 — An — 1) which ex- 
cludes n = 3. The case n = 2 is inconsistent with Eq. (|3*Tj) . unless e = +1, x < 0. 
For n = 5, 6, . . . the geodesies are timelike, for n = 4 it is a spacelike geodesic. 





a 


Eq. Pjl 


Eq. (jHU 




i(fc-n-l) 


4A 2 T 2 x fe + CT3X- 2 ss 0, fc < 


x 2 = 4A 2 ^x k + Cx~' 2 


|As6| 


-2 


1/x ~ for n = 2 


x 2 « 4A 2 /x 


|As7l 


-3 


1/x ~ for n = 4 


x 2 sa 12^ 2 /x 


|As8| 


-(n + 1) 


Cx" 2 ~ 


x 2 = -4^ 2 x"^ +1 ^ + Cx" 2 


|As9| 


-2n/(n - 1) 


Cx~ 2 ~ , n^3 


. 2 = 4av-i)» -^ir 1 +C x- 2 

(n+l)(n-3) 


Asio 


-A(n + 1) 


Cx" 2 ~ , n^3 


^ = + Cx- 2 



Table 3: Asymptotic solutions of the form u = Ax a as x — > 00 admit 0. All geodesies are 
null (e = 0), except in the last case lAslOl for which e / 0. 

The corresponding asymptotic forms of x(r) near x = 00 are summarized in table 0J 



8 





explicit solutions as x — > oo 


1 A , , -1 


x « 
x pa 


(2 - 

2VC(r-r )" 


4 (r - r ) 

i 

2 


2 

2 " fc for - 2 < fc < 
for k < -2 


IAs6l 


x w [3A(r-r )]3 


|As7l 


x pa 


3>/3A (r - T ) 1 


|As8| 


-x sa [(n + 3)A(r -r )]^s for C = 

i 

xpa 2VC(t-t ) 2 for C^O 


|As9| 


x sa 
x pa 


^ n — 6n+] 
2^(r-r ) 


_ 3) (r to) 
i 

2 


2(1— n) 
n 2 -6n + l 

for n > 3 
for n = 2, C / 


Asio 





Table 4: Explicit asymptotic geodesies as x(t) — ► oo. 

Let us note that there exist some special non-trivial exact solutions. For example, the 
particular case n = 2 of I As8l implies T3 = 0, see (|33[1 . and the normalization condition (|32j) is 
thus exactly satisfied. The corresponding null geodesies are given by 

/ x 3yl 2 \ 6j4 4 / 2t4 2 \ 

VCx2-4^ x ^_ + _j + __i n ^Cx 2 -4^2 x + x VC-- 7 =J =t-t , (37) 

where C > 4j4 2 /x > 0. For large x this approaches the asymptotic form presented in table 0J 
Similarly, the case n = 5 of I As 101 gives timelike geodesies (e = —1) with C > 0, 

x = ^/(t-t ) 2 ±2VC(t-t ) . (38) 

This resembles the asymptotic solution ()29(l but there is a different sign of the term (r — To) 2 . 

It follows from tables OH that the assumption (|30|) admits a large class of asymptotic null 
geodesies on the hypersurface y = 0, in particular in the region x — > 00 where the family is 
parametrized by four constants of integration, namely k, A,C,tq, and uq. The parameter a 
is in general negative (except for the possibility lAslj) so that u —* 00 as x —* 0, and -u — ► 
in the region x — ► 00. Therefore, u — > uq as x — > 00 which means that the geodesies become 
asymptotically tangent to the corresponding wave surface. Numerical simulations also indicate 
that near x — > the possibility IAs3l corresponding to a = — 1 (i.e. u ~ 1/x) is preferred. 

3.3 Exact "power-law" geodesies 

We have already mentioned that there exist three explicit families A n) , B n \ C <n) of exact solutions 
to the geodesic equations ([23 )1 - (|25 |) with G given by (|22|). These are particular subclasses C = 
ofEHEHilEIini respectively, cf. lf33 j) -(j3l) |) and table El In view of Eq. (j31[) . the coordinate 
x must be some power of the affine parameter r. It then follows from Eqs. (|3U|) and Q25|) that 
also the coordinates u and u must have the same behaviour, i.e. 

x = a{r- tq) p , u = (3(t - T ) q + u , v = 7(7- - T ) r . (39) 
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The constants a, (3,j,p,q,r,T ,u , related to previous parameters by 

a = (q-l)/p, A = qp/a a , (40) 

are restricted by the above equations. A straightforward calculation shows that there are three 
distinct possibilities, according to the value of the parameter B = f}y. 

A™' These solutions which exist for all powers n are given by B = —1, but a n_1 < so that this 
only apply to even n (implying x < 0) . The geodesies A n are null and correspond to IAs81 
characterized by a = — (n + 1), with the choice C = 0. 

B (n) The solutions B (n) arise when B = — 2/(n — l) 2 and correspond to the C = subcase of lAs91 
given by a = —2n/(n — 1). These geodesies are also null in the region x < for n = 2, or 
in the region x > if n > 3. 

C (n> The solutions C (n) for B = -±(n - 3)/(n - 1) in the same domains of x as B are spacelike 
(e = 1) for n = 2 or 4, and timelike (e = —1) for n > 4. This is the possibility lAslOl 
characterized by a = — (n + l)/2, for C = 0. The parameter 7 is now uniquely given as 
27 2 = (n — 3)a n ~ 1 . For n = 2 this implies a < 0, i.e. x < 0. 





B 


e 


V 


g = — r 


a"- 1 


/3 




-1 





2 


n — 1 


7 2 


1 


n + 3 


n + 3 


(n - l) 2 


7 


B ( '° 


- 2 





2(1 - n) 


(n- l) 2 


J_( n _ 3 )( n+ 1) 


-2 


(n - l) 2 


n 2 — 6n + 1 


n 2 — 6n + 1 


7(n- l) 2 


C (n) 


-(n-3) 


±1 


1 


n — 1 


n-l 

/ -4e y 


- ln-3 


2(n - 1) 


2 


\ n 2 - An - 1 y 


27 n-l 



Table 5: Three classes of exact solutions in the hypersurface y = 0. The families A , B of null 
geodesies are parameterised by an arbitrary constant 7. 

Finally, we also attempted to find more general geodesies such that y 7^ 0. Assuming 

x = a(r-r ) p , u = f3(r - T ) q + u , v = 7(r - r ) r , y = 5{t - r ) s , (41) 

we found that there are no such geodesies when n > 2. However, for n = 2 (i.e. = x 2 — y 2 ) 
there exists a class of exact null geodesies 

3 0/ \ 2 / \I 1/ \-I 3\/5 n. 2 

s =o7 (i"- T o) 5 , t; = 7(T-r )5 , n = (r - r ) 5 + u , y = ±-x~7 ( r ~ T o) 5 • (42) 

2 7Z 

4 Parallelly transported frames 

After analyzing privileged families of geodesies in the Kundt spacetimes we proceed to the 
study of singularities at x = and x = 00. Their character is not obvious: they cannot be scalar 
curvature singularities [28,29]. However, this leaves open the possibility that these are either 
non-scalar curvature singularities or quasi-regular singularities, according to the classification 
scheme introduced in [36]. To resolve this problem it is necessary to find orthonormal tetrads 
which are parallelly transported along geodesies approaching x = and x = 00. Projections of 
the curvature tensor onto these frames will then elucidate the character of the singularities. 
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4.1 Frames parallelly transported along any timelike geodesic 



A natural "interpretation" orthonormal frame {em(r)} for the Kundt spacetime ®, 

e (o) = « M = (v,x,y,u) , 

efij = -(i/(2x 2 t t ), 1,0,0) , (43) 

e f 2 ) = -(y/(2x 2 t t ), 0,1,0) , 

ef 3) = + (l/(2x 2 u), 0,0,0) , 

was introduced and employed in [31]. Mutually perpendicular spatial unit vectors em, em, er 3 \ 
are orthogonal to the four-velocity u of a timelike geodesic observer. Moreover, the vector e( 3 ) 
is the projection of the principal null direction k. Equation of geodesic deviation expressed in 
the tetrad (|43[). see Q61|l below, exhibits the transverse character of gravitational waves with e( 3 ) 
representing their direction of propagation (see [31] for more details). 

However, in general the tetrad (|43|) is not parallelly transported. Nevertheless, orthonormal 
frame {ell(j)(r)} transported parallelly along any geodesic must exist, and at each event it has 
to be related to the frame (|43jl by a 3-dimensional rotation of its spatial part. We may naturally 
assume that such relation is parametrized by the Euler angles <p, tp, 



' eli(i) ^ 




/ COS if 


e\ 2) 




sin ip 


V ell (3) J 




K 








/ COS l/j 


— sin ifi \ 


f e (i) ^ 


sin^ 


cos'i/' 


e (2 ) 


V 


1/ 


I e (3) / 



(44) 



Now, we impose the conditions for the parallel transport of the vectors ell(j)(r). After straight- 
forward but somewhat lengthy calculation, using (JHJ) and Eqs. (fTK|) - (|18|) . we obtain the following 
system of differential equations for the Euler angles, 



x sin 1? ip = — cos ifi , 

x sin $ tp = —y sin t9 + cos $ cos ip . 



(45) 



For a given timelike geodesic, characterized by specific functions x(t), y(r), we thus obtain a spe- 
cific rotation described by (p(r), $(t), ip(r) which define the parallelly transported orthonormal 
frame {ell(j)(r)} rotating with respect to (|l3*)) . 

Interestingly, the situation simplifies considerably for geodesies restricted to y = 0. The 
vector e%\ = (0,0,-1,0) is orthogonal to this hypersurface and it is parallelly transported. It 
only remains to find a suitable rotation of the vectors e^, e^, which may be obtained from 
(I45j) for a particular choice = vr/2, tp = 7T, implying 

(P = x-\t). (46) 

The parallelly transported orthonormal frame is thus given by 

e"(i) = — cos ip e (1) +sin^ e (3) , 

e ll( 2 ) = - sin if e (1) - cos ip e (3) , (47) 

e "(3) = ~ e (2) > ell (0) = u > 
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or explicitly 

e " (i) = cos V / {2x 2 u), 1,0, O^j — sin ip (v — l/(2x 2 ti), x, 0, uj , 

e^ 2) = simp (x/(2x 2 u), 1,0, o) + cosv? - l/(2x 2 ti), x, 0, , (48) 

ellf 3) = (0,0,1,0) , ellf 0) = (t>,x,0,u) . 

The vectors en) , e(3) of (|43|) — the latter indicating the direction of propagation of gravitational 
waves — rotate with respect to parallelly transported frames along any timelike geodesic, where 
the angle <p is obtained by integration of Eq. (|4(j|). Notice that <p — > oo as x — > 0. 

4.2 Frames parallelly transported along any null geodesic 

In order to find parallel frames along null geodesies it is more convenient to employ the formalism 
of null complex tetrads {k, 1, m, m} such that k • 1 = —1, m • rh = 1, see, e.g., [4]. We may start 
with the natural tetrad (j3J) and perform its boost followed by a null rotation with fixed 1, 

k' = Bk + Km + Km + KKB^ 1 1 , 

1' = B" 1 1 , (49) 
m' = m + KB x l , 

where the parameters are given by 

K = C, B = v-F/(2Q 2 )u , 2&-2Q 2 uv + Fii 2 = , (50) 

cf. 0. We arrive at the null tetrad 

k' = v d v + C + C + u d u , 

1' = FQ- 2 {2Q 2 v - Fii)- 1 d v + 2{2Q 2 v - Fii)' 1 d u , (51) 
m' = FQ~ 2 ((2Q 2 v - Fii)" 1 d v + + 2C{2Q 2 v - Fit)- 1 d u , 

for which the null vector k' is obviously tangent to the corresponding null geodesic and it is, of 
course, parallelly transported along this. We wish to find a tetrad such that the remaining null 
vectors are also parallel. This is achieved by performing additional null rotation — this time 
with fixed k' — and a spatial rotation in the transverse plane, 

k = k' , 

HI = l'lLm' + Lm'l LL k' , (52) 
mil = exp(i0) (m' + Lk') . 

We have to choose the complex function L and the real function of the affine parameter 
such that 111 , mil and mil are parallelly transported. Calculation using (|50j) and the Christoffel 
symbols for Q leads to the system of differential equations 



Q 2CC 



+ — WC^c - C^,c) - (C - OH) , 



L = + i[Q(CH ( - C# <) - (C - C)H}} L-^~ —dQH,c - H) . (53) 

^ (J OAA ' s > U Iff ' s 



Q 2CC J Q 2CC 
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In real coordinates (0) of the metric (JSJ these conditions take the form 

y iu 2 

<P = To ; -2 [ x ( ±G ,y ~ V G ,x) + V G \ 

x x + y 



2\/2u 



2 



L x — L 2 <p = -TY——2 {xG x — G) j=— , (54) 



- , - 'J 2 

L 2 + Li 6 = — = ^ xGm , 

x 2 + y 



where L\,L 2 denote real and imaginary parts of L, L = L\ + i-^2- These equations are more 
complicated than the corresponding expressions 1)4 5 jl for parallel frames along timelike geodesies. 
In particular, they explicitly depend on the structural function G. 

However, expressions (|54[) again simplify for null geodesies on the hypersurface y = in the 
Kundt spacetimes (j22j) . The simplest parallelly transported frame is then obtained by setting 

cf) = 0, L 2 = 0, U = 2y/2{n- l)x n x~ 2 u 2 - (v^x) -1 . (55) 

In real coordinates x, y, the frame is thus explicitly given by 

fell** = (v,x,Q,u) , 

= {l\v + 2iix~ 2 {v 2 + G/x)(l + \Z2Ltx), L\x + \f2L X} 

0,L 2 1 u + 2ux- 2 (l + V2L 1 ±)^j , (56) 
m f = (t, x v + y/2 iix' 1 ^ 2 + G/x),Lix + 1/V2, 0, L x u + V2ux 
mf = (0,0,l/\/2,0) . 



where mil = m| + im||. For particular null geodesies of the form (|3Uj) given by equations Q31|). 
?2|) with e = 0, C = we obtain 



Li = ^ 5 (n + 2a + 2)x- 1 (r) . (57) 
This resembles analogous expression (|46(l for timelike geodesies. 

5 Frame components of curvature tensor and geodesic deviation 

Now we evaluate explicitly the components of the Riemann tensor in the above tetrads which 
are parallelly transported along timelike or null geodesies in vacuum Kundt spacetimes. 

5.1 Curvature components evaluated along timelike geodesies 

Considering we first project the Riemann tensor onto the "interpretation" frame (|43|). 



#(0)(1)(0)(2) 



%(1)(0)(1) 


= %)(1)(1)(3) 


= #(1)(3)(1)(3) = 


A+ . 


#(0)(2)(0)(2) 


= #(0)(2)(2)(3) 


= #(2)(3)(2)(3) = 


-A+ 


= #(0)(1)(2)(3) 


= #(0)(2)(1)(3) 


= #(1)(3)(2)(3) = 


-A x 



( 58 ) 

) = #(0)(1)(2)(3) = #(0)(2)(1)(3) = #(1)(3)(2)(3) = ~A X , 

where the amplitudes are 

A+ = -2xii 2 G >xx , A* = 2xu 2 G, xy , (59) 
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see also [31]. A straightforward algebraic calculation using the transformation (|44j) leads to the 
following non-vanishing components of the curvature tensor in the orthonormal frame {ell(j)(r)} 
parallelly transported along a timelike geodesic, 



-R 
-R 
-R 

R 
-R 
-R 

R 
-R 

R 
-R 

R 

R 



(2)(3)(2)(3) 
(0)(1)(2)(3) 
1)(3)(1)(3) 
(0)(2)(1)(3) 
1)(2)(2)(3) 
(0)(3)(2)(3) 
(0)(3)(1)(3) 
1)(2)(1)(3) 
1)(3)(2)(3) 
(0)(2)(2)(3) 
1)(2)(1)(2) 
(0)(3)(1)(2) 



A' + (cos 2 p — sin 2 93 cos 2 i?) + A' x cos $ sin 2<£> = ^L)(i)(o)(i) 
A' x (cos 2 tp — sin 2 ip cos 2 $) — A' + cos $ sin 2<p , 

(sin 2 <^9 — cos 2 ip cos 2 "!?) — A' x cos 1 !? sin 2ip = ^L)(2)(o)(2) 
A' x (sin 2 ip — cos 2 (p cos 2 "!?) + A' + cos $ sin 2(p , 
A' + i sin 2$ sin p — A' x sin -d cos <p = R] 



2 ' 
\i 1 



L (0)(l)(0)(3) 

-R\ 



L (0)(l)(l)(2) 
(0)(2)(1)(2) 
(0)(2)(0)(3) 



A' x 2 sin 2-d sin 99 + A' + sin $ cos </? 
A' + sin $ sin 99 — A' x \ sin 2$ cos p 
A' x sin $ sin ip + ^ sin 2d cos <£> 
•4+|(l + cos 2 "9) sin 2</3 — .A x cos $ cos 2y? = ^L)(i)(o)(2) 
.4 X ^(1 + cos 2, d) sin 2c/? + .A' + cos $ cos 2<^ = ^L)(i)(i)(3) 
A' + sin 2 tf = 
A' x sin 2 i? , 



(60) 



t (0)(3)(0)(3) 



where the Euler angles satisfy (J32J), A' + = A+ cos 2^+^4 x sin 2ip, A' x = —A+ sin 2^+-A x cos 2ip. 

The above components of the Riemann tensor can be employed to study geodesic deviation. 
Let us consider a displacement vector Z^(r) between two nearby particles moving along timelike 
geodesies. We project Z^ onto {e^(r)}, Z® = effz 1 *, and denote its absolute derivatives as 
= e^(DZ^/dr), = e]L (D 2 Z^/d 2 r). The equation of geodesic deviation with respect to 



(P|) is Z® = -R 



(0 

(o)CO(o) 



Z®, i.e. using 



-A + ZU+A X Z& , 
A+ZW+A X ZW, 



(61) 



. 



see [31]. The structure of the right-hand side exhibits the transverse character of gravitational 
waves with two polarization modes (represented by the amplitudes A+,A X ) which propagate in 
the spatial direction of e( 3 ) . However, the left-hand side is complicated because the frame (|43[) 
is not parallelly propagated. Nevertheless, using the relation we can express 111 

terms of the derivatives with respect to r. The equation (|61j) thus takes the form 



d 2 ZW 2y dZW 2 dZ® 
dr 2 x dr x dr 



d 2 Z( 2 ) 2ydZW 



( A+ -\-t;)zU 



+ 



(Ay 



+ 



X 

yx 

x 2 



rZ (3) 



dr 2 



x dr 



Ay 



X x' 



(62) 



+ (A + + ^)ZM-A;ZW 



II 
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d 2 Z^ 2dZ« x m y m 1 m 
dr z x dr ar x z x z 

In the Kundt spacetimes given by G< n > there are timelike geodesies restricted to the privileged 
hypersurface y = 0, see l|22|). In such a case we obtain «4 X = 0, 1? = 7r/2, = it, and 99 solving 
(|46|). The only components (j6()fl of the Riemann tensor with respect to the frame (|47fl are thus 

_ -^Jwairarai = -4+ cos2 V = Rl 



L (2)(3)(2)(3) "(0)(1)(0)(1) ' 

_jR (l)(3)(l)(3) = -4+ sm V 3 = - R (0)(2)(0)(2) ' 

_jR (0)(3)(2)(3) = "4+ COS ^ = _jR (0)(l)(l)(2) ' 

- R (0)(3)(l)(3) = Sill< ^ = _jR (0)(2)(l)(2) ' ( 63 ) 

jR (l)(3)(2)(3) = -4+ sin (/? COS (/? = i?( )(l)(0)(2) ' 

i? (l)(2)(l)(2) = A + = - R (0)(3)(0)(3) ■ 

The equations of geodesic deviation (ff)2"|) also considerably simplify since y = = .4. x • 
5.2 Curvature components evaluated along null geodesies 

We also calculate components of the Riemann tensor in the parallelly transported null tetrad 
(|52|). For vacuum spacetimes the curvature tensor coincides with the Weyl tensor. It is conve- 
nient to use standard transformations of the NP coefficients tyj, j = 0, 1, 2, 3, 4, see e.g. [4]. We 
start with the expression Q with respect to the tetrad |2J), and perform boost B, null rotation 
K, followed by null rotation L and spatial rotation (j), where these parameters are given by ()5fl|) 
and (|53[). respectively. We obtain the components, 



= e (2-j)i0 R 4-j ^ + RL y y ^ ^ 



where, using (J5U|) 



2 



2 



2(CC) 

8xu 2 Af + i4. x 

~~ 7 • 2 1 r 2\l2 ' — ",z:z + iG xy ) — 4 — - . (65) 

(ar + y z ) z (x z + y z ) z 

Comparing (|64[) with Q6UJI for timelike geodesies (for which only A+ and A x are sufficient 
to describe the behaviour of gravitational filed), investigation of the curvature tensor along null 
geodesies is more complicated. One has to study not only ^ but also various powers of K and 
(1 + KL). However, there is a simplification for null geodesic y = in vacuum Kundt spacetimes 
with G^; we obtain * = 4A + /x A , (p = 0, K = x/y/2, and L = L\ given by lfHK|) . so that 



*]• = (Li + V2/iyA+ , 3= 0,1,2,3,4. (66) 



6 Character of the singularities in the Kundt spacetimes 

Using the above results we will now discuss the behaviour of the curvature tensor components 
with respect to frames parallelly propagated along timelike or null geodesies. This will elucidate 
the character on the singularities, in particular of the envelope singularity x = 0. 
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It has been shown that such components of the Riemann tensor are proportional to the 
amplitudes A+, A x , see expressions l|50*|l or for timelike or null geodesies, respectively. 
We calculate these amplitudes (|59*|) for spacetimes given by the structural function (|2*Tj) with 
c = 1 on the hypersurface y = as A+ = — 2n(n — 1) x n_1, u 2 , A x = 0. The behaviour of A+ as 
x — > or x — > oo depends on the mutual relation between u and the powers of x. This again 
justifies our assumption PU|) which enabled us in section El to derive particular classes of exact 
or asymptotic geodesies, summarized in tables HOD For all these we obtain 

A + = -2n(n - l)A 2 x 2 ^- 1 . (67) 

It is now easy to express explicitly the behaviour of this function for exact geodesies denoted 
above as A"\ B (n) , C (n) . Recall that only the geodesies C <n) are timelike if n > 4, and the behaviour 
of A+ is thus sufficient. All other geodesies are null and for the discussion of (|fi6[) it is necessary 
to evaluate also the function L\ + \/2/x. In table we present a summary of the powers of 
(r — tq) corresponding to each quantity. 



Class 


e 


X 


A+ 


L x + y/2/x 


*J = (Li + V2/x) 3 A + 


A n) 





2 


-2 


n + l 


(i-2)n + (j-6) 


n + 3 


n + 3 


n + 3 


B (n) 





2(1 -n) 


-2 


n 2 — An — 1 


C?-2)n 2 -4(j-3)n-(j + 2) 


n 2 — 6n + 1 


n 2 — 6n + 1 


n 2 — 6n + 1 


c (n) 


-1 if 

n > 4 


1 


-2 







Table 6: The powers of (r — tq) for various quantities which enter the components of the cur- 
vature tensor in frames parallelly propagated along classes of exact geodesies y = 0. 

Interestingly, although the power-law dependence of x ~ (r — tq) p on the affine parameter r 
varies for different classes A n \ B ( "', and C <n) of geodesies, the dependence of the curvature tensor 
amplitude A+ is the same, namely 

A+ ~ (r - T r 2 . (68) 

This is the consequence of the fact that 2a + n — 1 = 2(g — l)/p + (n— 1), cf. (|40[) . and also 
q = —(n — l)p/2, which can be seen from table 03 so that 2a + n — 1 = —2/ p. Similarly, the 
relation (|57|) gives L\ ~ ~ (r — to) _p , implying Li ~ (r — to) 1 ^^. Since cc ~ (r — to) p_1 , 
both Li and 1/x exhibit the same behaviour, 

L 1 + V2/x~(T-T ) 1 - p . (69) 

Let us now discuss specific behaviour of the curvature tensor in some more detail. 

The only class which admits timelike geodesies is C when n > 4. This is characterized by 
a = — (n + l)/2 and p = 1, i.e. x ~ (r — to). The amplitude (|55|) thus behaves as 

lim A+ ~ lim (r — to)~ 2 = oo , 

X— >0 T-»T0 

lim ^4+ ~ lim (r - r )~ 2 = . (70) 
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Obviously, all non-trivial components (|63|) of the curvature tensor with respect to frames par- 
allelly transported along these timelike geodesies approaching x = diverge. The envelope 
singularity of the radiative Kundt spacetimes is thus a non-scalar curvature singularity which 
is reached in a finite value of the proper time. On the other hand, the spacetimes become flat 
in the region x — > oo for these timelike geodesies. 

Analogous discussion concerning null geodesies is more complicated since the components 
(|66|) depend on A+ and also on powers of the function (L\ + \f2jx). The class of exact null 
geodesies A (n) is characterized by a = — (n + 1) and p = 2/ (re + 3), so that x ~ (r — tq) 2 ^ u+3 \ 
A+ ~ x~( n+3 ) ~ (r - t )- 2 . The relation (MJ) reads (L x + y/2/i) ~ (r - To )(«+i)/(n+3)_ Thc 
frame components of the curvature tensor are thus 

#J ~ ( T _ To )[(i-2)™+(i-6)]/(n+3) _ ( 71 ) 

The lowest power appears for j = 0, namely = A+ ~ (r — to) -2 , whereas the highest one 
for j = 4, ^| ~ (r — To) 2 ( n ~ 1 - ) /( n+3 ). Obviously, there is always a curvature singularity at x = 
since the frame component of the curvature tensor ^ diverges, 

lim M ~ lim (r - r ) -2 = oo . (72) 
Other components also diverge at x = if the corresponding power in (|7T|) is negative, i.e. for 

i<2 + — L (73) 
re + 1 



Thus, ^q, always diverge, ^ diverges at x = for spacetime with re = 2, whereas \I/ 4 

always approaches zero there. For null geodesies of this class A (n) approaching the region x — > 00 
the situation is completely the opposite: Wq, Vl/j, ^ always vanish asymptotically, whereas 

Jim # 4 ~ lim (t - To )2(«-i)/(«+3) = oo _ (74) 



Behaviour of the curvature tensor along null geodesies B (n) with a = —2n/(n — 1), p = 
2(1 — n)/{n 2 — 6re+l), depends on the parameter n of the specific function (|21|) of the spacetime. 
In particular, x ~ (r — ro) p so that 

for n<6: x — > 44> r ^ To , x — > 00 44> t ^ 00 , 

for n > 6 : x ^ 44> r ^ 00 , x ^ 00 44> t ^ to . (75) 

The behaviour is thus more complicated. The amplitude A+ = ^ is proportional to (r — to) -2 , 
i.e. singular as r — * tq, other components of the curvature tensor are 

\$\\ ^ ( T _ To )[(j-2)n 2 -4(j-3)n-(j+2)]/(n 2 -6n+l) _ ^ 

The frame component diverges as r —* ro if the corresponding power of (|76j). which is equal to 
j(l — p) — 2, is negative, i.e. (considering p < 1) when 

2 4(re - 1) 

i<! = 2 ~ 2 A \ - ( 77 ) 

1 — p re — 4re — 1 

The behaviour for r — > 00 is exactly the complementary one. In particular, for spacetimes with 
n = 2 the components *$>q, *$>\, ^ diverge at x = whereas Vl/3, ^ 4 approach zero there. For 
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n = 3 the components Wq, w[, W^^i diverge at x = while is regular, for n = 4 all the 
components diverge. In spacetimes n = 5 with p = 2 the condition reads j > —2 so that all 
components ^/j diverge as x — > 0, but all approach zero as x — > do. For n = 6 the condition 
(|77|) is satisfied only for j = so that, considering (f?5|) . VP^, ^35 ^4 diverge at x = 0, and 
diverges as x — > 00. For spacetimes with very large parameter n the condition is satisfied when 
j = or j = I so that ^3, ^4 diverge at x = 0, and ^g, ^1 diverge as x — > 00. 

The above described behaviour of the curvature tensor near x = and x = 00 has been 
obtained using explicit privileged classes A™ , B (n) , and C (n) of exact geodesies on the hypersurface 
y = 0. This can be confirmed by (the only) "power-law" geodesies of the type (|41j) which admit 
y 7^ 0. These null geodesies in spacetimes with n = 2 are explicitly given by (|42|) . There is 
x ~ (r — t ) 2 / 5 ~ y, u ~ (t — r ) -1 / 5 so that, according to (JH§j), *4x = 0, .A+ ~ xit 2 ~ (r — r )~ 2 , 
and K ~ (r — to)" 3 / 5 , (1 + XL) = const. Consequently, from (|64|) we obtain 

^~(r-ro) 3j / 5 - 2 , (78) 

which implies a divergence of the components ^q, ^31 a ^ # = 0. On the other hand, 

^4 ~ (t — r o) 2//5 ~ which indicates a singularity at x —* 00. 

Finally, we investigate the behaviour of the curvature tensor along more general classes of 
asymptotic geodesies approaching the regions x = or x = 00. These have been described in pre- 
vious section EH! The explicit class of null geodesies lAsll is characterized by x ~ (r — tq) 2 ^ 2 ^^ 
where k = 2a + n + 1 , see tables ^ and |2j In view of (|67|) we obtain 

A + ~x k ~ 2 ~{t-t )- 2 , (79) 

so that the component = ^+ diverges at x = as r —* tq. Other components are given by 
~ (r — To)"i~ 2 ) fc + 4 ]/( fe-2 ) . The same applies to asymptotic geodesies IAs2l which is a particular 
subcase of lAsll for k = n — 3. The classes IAs3l and lAs4l exhibit somewhat different dependence. 
It turns out that ^ ~ (r — tq)^ -3 "^" 2 and ~ (r — ro)( n-7+ ^' 2 3 respectively. Therefore, as 
x — > with r — > ro, only the component \& ' diverges for n = 2 in case I As3l and for n = 6 in case 
IAs4l Other components of the curvature tensor are regular at x = 0. 

Concerning the complementary asymptotics x — ► 00, for the case IAs5l the functional depen- 
dence of x is the same as in the case lAsll fcf. tableHJ so that the behaviour of the amplitude is also 
given by ()79|) , With r — > 00 the component \I/q vanishes, but \&'! diverge when j > 2 — 4/k > 2. 
Since lA"s6l and lA"s71 are particular cases of IAs5l for k = —1, all the components vanish asym- 
potically along these special geodesies. Finally, there are three subclasses IAs8| [As91 lAsTUl For 
C = 0, x — > 00, they asymptotically approach the exact geodesies A (n) , B (n) , C (n) which we have 
discussed in detail at the beginning of this section. For C / there are particular null geodesies 
of the form x ~ \/t - to, for which *!j ~ (r - T ) (j ~ n ~ 3)/2 in case lAs81 and *|J ~ (r - r ) (j_7 ^ 2 
in case IAs9l In both these special cases, all components vanish as x — > 00, i.e. as r — > cxd. 

To summarise our somewhat lengthy discussion: typical components (like A+ = ^q) of the 
Riemann tensor with respect to parallelly propagated frames diverge as x — > (which is reached 
for r — ► To) froi/i /or timelike and null geodesies, see relations (J70J), (|71 j) -([73 j) . (|76|). (|78|1. and 
(|79|l. Concerning the behaviour in the complementary region x — > cxd of the Kundt spacetimes 
there are also divergences in other components of the curvature tensor (typically as r — > cxd. 
However, this is not so along all null geodesies and, in particular, along timelike geodesies of the 
class C (n> for which A.+ — > 0, see ()7U|) . Further analysis is thus required to fully understand the 
global structure of null and timelike infinities of these spacetimes. 
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7 Conclusions 



In the present work we have studied vacuum solutions to Einstein's equations which represent 
exact radiative spacetimes of the Kundt type. Our main objective was to investigate the geodesic 
motion and to elucidate the character of the singularities contained. 

We found, in section 01 particular families of exact geodesies denoted as A (n> , B <n) , C <n) (cf. 
table 121, and also classes of approximate solutions to the geodesic equations on the hypersurface 
y = as x — > ()Asll to fAs4|) or x — > oo (|As5l to lAslOjl . see tables El and |U These asymptotic 
families are quite large as they contain several independent parameters, e.g., k, A, C, tq, uq. 

In section |1] we presented tetrads that are parallelly transported along any timelike or null 
geodesic. Somewhat surprisingly, these could be given in a closed explicit form. In particu- 
lar, we demonstrated that the spatial direction e(3) of propagation of the Kundt gravitational 
wave always rotates with respect to parallelly propagated frames. This is similar to analogous 
behaviour which was discovered previously in the Siklos spacetimes with A < [37]. It can be 
explained by the fact that both Kundt and Siklos waves belong to a larger class of generalized 
Kundt spacetimes for which the spin coefficient r is non- vanishing [19], see also [5]. 

Subsequently, in section [3 we calculated projections of the curvature tensor onto the frames 
that are parallelly propagated along arbitrary timelike or null geodesic. Not only could these 
components be used for studies of geodesic deviation, see equation (|62|) . but they are also crucial 
for discussion of the character of singularities. 

This is done in the final section H3 From the frame components of the curvature tensor 
we explicitly demonstrated that there is a physical singularity at x = which is reached in 
a finite value of the affine parameter tq both for timelike and null geodesies. Although all 
scalar curvature invariants of these Kundt spacetimes identically vanish, the expanding envelope 
of the rotated wave surfaces u = const, (see figure that is localized on x = is a (non- 
scalar) curvature singularity, according to the classification scheme introduced in [36]. It is 
not a "mild" quasi-regular singularity because some of the Riemann tensor components in the 
parallelly propagated orthonormal frames diverge as x — > 0. It is thus not possible to extend 
the Kundt solution "smoothly" into the inner part of the expanding wave-front envelope x = 
by a simple flat Minkowski region. 

The complementary region x — > oo of the Kundt spacetimes corresponds to conformal in- 
finity because it is reached by geodesies only as r —* oo. There are also divergences in some 
components of the curvature tensor, even if not along all geodesies. A more detailed analysis 
is still necessary to understand the global structure of these null and timelike infinities of the 
Kundt wave spacetimes. 

Finally, let us observe that many of the above results obtained for type A?" vacuum Kundt 
waves remain valid also for conformally flat spacetimes © containing pure radiation. In partic- 
ular, on the hypersurface y = the metric function Gq = x 2 + y 2 (including its first derivatives) 
is effectively identical to the function (|21|) with n = 2, i.e. 

G (2) 

= x 2 - y 2 . We thus obtain 

the same forms of the geodesies and parallelly propagated frames. The tetrad components of 
the curvature tensor (|63j) are also similar, but expressions (|64|) now vanish identically because 
the spacetimes are conformally flat (\E r = 0). However, projecting directly ()1U|) onto ()56|) we 
obtain, for example, R a ^i,m\ k^m^k^ = —2xu 2 = — 2A 2 x 2a+1 . For the class A (2> of exact 
null geodesies one has a = —3, x ~ (r — To) 2 / 5 , for B 2 there is a = —4, x ~ (r — To) 2 / 7 . In both 
cases the above frame component of the Riemann tensor is proportional to (t — To) -2 and thus 
diverges as x — > with r — > tq. Obviously, the envelope singularity x = in conformally flat 
Kundt spacetimes with pure radiation is again a non-scalar curvature singularity. 
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